Is the term "type-1.5 superconductivity" warranted by Ginzburg-Landau theory? 



O 

(N 

C 

m 



c ■ 
o ■ 
o 

u ■ 

^ : 

c3 ■ 



I 

c 

o 
o 



(N 
> 

00 

in 
p 

00 

o 

O 



X 



V. G. Kogan, J. Schmalian 
Ames Laboratory and Department of Physics & Astronomy, Iowa State University, Ames, I A 50011 

(Dated: Aug.2, 2010) 

It is shown that within the Ginzburg-Landau (GL) approximation the order parameters Ai (r,T) 
and A2 (r,T) in two-band superconductors vary on the same length scale, the difference in zero-T 
coherence lengths £01/ ~ fvupv /A„(0), v = 1, 2 notwithstanding. This amounts to a single physical 
GL parameter k and the classic GL dichotomy: k < l/y/2 for type-I and k > l/y/2 for type-II. 

PACS numbers: 74.25.Nf,74.20.Rp,74.20.Mn 



I. INTRODUCTION 

The physics of all superconductors near their critical 
temperature, T c , is based on the GL theoryj 1 - This in- 
cludes multi-band superconductors with distinct sheets 
of the Fermi surface. A number of recent papers deal 
with two-band materials with coefficients of the GL free 
energy (for the field- free state), 



F= (a»K + buA 4 j2)-2 1 A 1 A 2 . 



(1) 



introduced phenomenologically, see, e.g., Ref.[2|. Choos- 
ing these coefficients in various ways, one could arrive to a 
number of the choice dependent conclusions. 3,4 However, 
the coefficients can be derived from microscopic theory; 
they are certain functions of the microscopic coupling 
constants responsible for superconductivity and of tem- 
perature T. This has been done time ago by Tilled and 
later by Zhitomirsky and Dac 6 who have shown, within a 
weak-coupling model, that the coefficients a v do not have 
the familiar GL form a(T — T c ). Instead, they acquire 
a constant part, const + a(T — T c ), which is intimately 
related to the constant 7 of the mixed Josephson-type 
term to ensure A u oc \JT C — T near T c . 

We argue in this work that the ratio of the order pa- 
rameters is T independent in the GL domain, 



Ai(r,T)/A a (r,T) = const, 



(2) 



with the constant depending on interactions responsible 
for superconductivity. The one-dimensional version of 
Eq. © has first been obtained while solving the GL prob- 
lem of the interface energy between superconducting and 
normal phases relevant for the distinction between type- 
I and type-II two-band superconductors^ For the strong 
intraband scattering (the dirty limit), the result ^ has 
been obtained by Koshelev and Golubov provided the 
interband scattering could be disregarded 4 Here, we es- 
tablish this result for any problem in the GL domain. 

We show that the equations for Ai(r, T) and A 2 (r, T) 
are reduced to one independent GL equation. In other 
words, there is a single complex order parameter describ- 
ing the two-band superconductor in the GL domain and, 
as a consequence, a single length scale £ for spatial vari- 
ation of both Ai(r,T) and A 2 (r,T). 



Our results, along with the earlier critique^ and a com- 
prehensive review by Brandt and Das, 10 question validity 
of publications discussing properties of MgB 2 within the 
GL framework where each band is attributed with its own 
coherence length and sometimes even with its own pene- 
tration depth, see, e.g., Ref.HTI and references therein. 

We stress that our claim that the gap functions 
A„(r,T) change on the same length scale relates ex- 
clusively to the temperature domain, however narrow 
it could be, where the GL theory is valid. Out of this 
domain and at low temperatures in particular, different 
length scales ~ TivfvI A„(0) may enter and result in prop- 
erties substantially different from those in the GL region. 
Still, as long as the GL energy functional is used, the as- 
sumption of two coherence lengths cannot be justified. 

Below, we discuss the phenomenologic two-band 
GL theory and later confirm our conclusions within a 
weak-coupling microscopic scheme. 



II. TWO-BAND GL IN FIELD 



The two-band GL functional reads: 
T=Uv{Y J U\A V ? + \\A V \ 

v=l,2 \ 



7 (a 1 a; + a 2 aj) 



8tt J 



(3) 



where II = V + 2iriA/(fro and the constant 7 along with 
the coefficients a, b, K will be given later. The GL equa- 
tions are minimum conditions for the functional (j3|). One 
obtains varying T with respect to A*: 



aiAi + feiAi|Ai| 2 - 7A2 - JTirfAi = 0, 



(4) 



a 2 A 2 + 6 2 A 2 | A 2 | 2 - 7A1 - /f 2 n 2 A 2 = . (5) 

We now recall that in the one-band GL equation, 

aA + 6A|A| 2 - i^n 2 A = 0, (6) 

all terms are of the same order (1 — T/T c ) 3 / 2 = t 3 / 2 
(A cx t 1 / 2 , a oc t, and n 2 cx £~ 2 cx t). This is not so for 
Eqs. ((4]), ([5]) because 7 is a constant and a„ may contain 
constant parts. 
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Having this in mind, we express A2 in terms of Ai 
from Eq. (@} and substitute the result in Eq. ([5]) keeping 
only terms up to the order r 3 / 2 : 

(oioa - 7 2 )Ai + (6io 2 + 6 2 a?/7 2 )Ai|Ai| 2 

-(011:2+02^1)^1 = 0. (7) 
Similarly, one obtains an equation for A 2 : 

(oio 2 - 7 2 )A 2 + (b 2 ai + 6ia^/7 2 )A 2 |A 2 | 2 

- ( ai K 2 + a 2 K 1 )U 2 A 2 = 0. (8) 

In zero field, one has A 2 oc (aia 2 — 7 2 ), so that at T c , 
aia 2 — 7 
parts 



J — 0, and therefore a v must contain constant 



(9) 



(10) 
(11) 



such that ai c a 2c = 7 2 . 

Eqs. (|7|) and © for A„ can now be written as: 

- arAi + ^ Ai| Ai| 2 - iOI 2 Ai = , 

-arA 2 +/3 2 A 2 |A 2 | 2 -iai 2 A 2 = 0, 

with 

a — ai + a 2 , jK" = ai c i^2 + a- 2c K\ , 

Pi = ha 2c + b 2 a\ c h 2 , /3 2 = b 2 a lc + haUr* . (12) 

We note that within the accuracy of the GL theory, up to 
0(t 3 / 2 ), the equations for Ai and A 2 are coupled only 
via the vector potential. 

In particular, in zero field we have 



A,, 



or/ 0„ , 



so that the ratio 



A 2 (T) _ h 
A? n (T) A 



(13) 



(14) 



*20\ x ) Pi 

comes out to be T independent in the GL domain. 

Furthermore, one easily checks that for any solution 
A x (r, T) of Eq. {TO]), Eq. (dU is satisfied by 

A 2 (r,T) = Ai(r,T)VA//32. (15) 

In particular, this implies that in equilibrium Ai (r,T) 
and A 2 (r,T) must have either the same phases or the 
phases differing by ttf^- It is found in Ref.@ that for 
small 7 the ratio A 2 /Ai changes away of T c ; we note, 
however, that this deviation is beyond the GL accuracy. 
Reliable results beyond GL can be obtained only within 
microscopic approaches like Gor'kov or Bogolyubov - de 
Gennes theories. 

Moreover, introducing the order parameters normal- 
ized on their zero-field values, 

Ai A 2 

A^M = A2^r) = *' (16) 

both Eq. (jlpp and ([TTj) are reduced to one: 

K 



*(1 - |*| 



-n 2 *. 



(17) 



Thus, the length scale of the space variation of both Ai 
and A 2 , the coherence length, is given by 



i 2 = Kjar . 



(18) 



III. MICROSCOPIC WEAK-COUPLING 
TWO-BAND MODEL NEAR T c . 

To establish connection of GL equations with the two- 
band microscopic theory we turn to a weak-coupling 
model for clean and isotropic materials (not because 
these restrictions are unavoidable, but rather due to the 
model simplicity). 

Perhaps, the simplest formally weak-coupling ap- 
proach is based on the Eilenberger quasiclassical formula- 
tion of the Gor'kov equations valid for general anisotropic 
order parameters and Fermi surfaces^ Eilenberger func- 
tions /, g for clean materials in zero-field obey the system: 



= Ap - huf , 



9 



1 



(19) 
(20) 



A(Jfe) = 2ttT7V(0)^ /V(k,k')f(k',uj)) . (21 



Here, k is the Fermi momentum; A is the order parame- 
ter that may depend on the position k at the Fermi sur- 
face. Further, N(0) is the total density of states (DoS) 
at the Fermi level per spin; the Matsubara frequencies 
are given by hw — ixT(2n + 1) with an integer n, and wd 
is the Debye frequency; (...) stands for averages over the 
Fermi surface. 

Consider a model material with the gap given by 



A(k) = Ai, 2 , k e Fi, 2 , 



(22) 



where Fi,F 2 are two sheets of the Fermi surface. The 
gaps are assumed constant at each band. Denoting DoS 
on the two parts as i\Ti 21 we have for a quantity X con- 
stant at each Fermi sheet: 

(X) = (XiiVj + X 2 N 2 )/N(0) = n x X x + n 2 X 2 , (23) 

where 711,2 = Ni, 2 /N(0); clearly, n\ + n 2 = 1. 
Equations ([Ti?]) and ([20)1 are easily solved: 



/„ = A ,,//?„, g v = hw/p u , Pi 



A, 



2, ,2 



h 2 uj 



(24) 



where v = 1,2 is the band index. The self-consistency 
equation (|2"Tj) takes the form: 



A„ = ^ JfyA^A 
[i=i,i 



2vrT 



(25) 



where A„ M = N(0)V Vfl are dimensionless effective inter- 
action constants. The notation commonly used in litera- 
ture, Xvfi = n^Xjj^, includes DoS'. We find our notation 
convenient since, being related to the coupling potential, 
our coupling matrix is symmetric: X v ^ — A M „. 

It is seen from the system (|25[) that Ai j2 turn zero at 
the same temperature T c unless Ai 2 = and equations 
(f2l)j) decouple, the property that has been noted in earlier 
work^iS As T — > T c , A 1>2 0, and (3 -> hw. The sum 
over ui in Eq. (|25[) is readily evaluated: 



S = J2 



2ttT 



In- 



2huj]j 



In. 



D 

1.762V 



(26) 
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7 = 0.577 is the Euler constant. This relation can also 
be written as 



1.76 T c = 2hw D e' s . 

The system (|25l) at T c is linear and homogeneous: 

Ai = 5(niAnAi + n 2 A i2 A 2 ) , 
A 2 = S*(niAi 2 Ai + n 2 A 22 A 2 ). 

The zero determinant gives S and, therefore, T c : 

S 2 nin 2 r] - S(n\\vL + n 2 \ 22 ) + 1 = 0, 
rj = AnA 22 — A 12 . 

The roots of this equation are: 



(27) 



(28) 

(29) 
(30) 



S : 



niAn + n 2 X 22 ± yj (niAn - n 2 X 22 ) 2 + 4mn 2 A 



2n 1 n 2 n 



(31) 



Various possibilities that arise depending on values of 
are discussed, e.g., in Refs. [l4lil8l Introducing T- 
independent quantities, 



Si = A 22 - n x nS , S 2 = An - n 2 r)S , 



we write Eq. (1291) as 



S\S 2 



\ 2 

A 12 ' 



(32) 



(33) 



the form useful for manipulations below. 

If X12 = 0, Eq. (f3~T|) provides two roots 1/niAn and 
l/n 2 A 22 . The smallest one gives T c , whereas the other 
corresponds to the temperature at which the second gap 
turns zero. We note that this situation is unlikely; it im- 
plies that the ever present Coulomb repulsion is exactly 
compensated by the effective interband attraction. 

Since the determinant of the system (f2"5| is zero, the 
two equations are equivalent and give at T c : 



A2 
Ai 



1 — niXuS 



(34) 



When the right-hand side is negative, A's are of opposite 
signs. Within the one-band BCS, the sign of A is a mat- 
ter of convenience; for two bands, Ai and A 2 may have 
equal or opposite signs »^ 

After simple algebra, Eq. (|3~4"f can be manipulated to 



A, 
Ax 



Si 

s 2 



(35) 



We thus obtain by comparing with Eq. (TT4")) or ([15)) the 
ratio of phenomenological coefficients in terms of micro- 
scopic couplings: /3i//3 2 = S\/S 2 . We have seen above 
that within the GL approximation this ratio remains the 
same at any T in the GL domain not only for a uni- 
form field-free state (or for 7 — > 00 as in Ref.l20() but for 
any situation with A's depending on coordinates in the 
presence of magnetic fields. 

We note that the proportionality of Ai and A 2 has also 
been shown to hold within microscopic weak-coupling 



theory in the dirty limit by Koshelev and Golubov^ It is 
also worth mentioning here that the above proof of this 
proportionality based on the GL approach is quite gen- 
eral and holds for any scattering, gap anisotropics etc. 

In the following we use the GL coefficients obtained in 
Refs.HI and[^. In our notation they read: 



7 



N(0) 
V 

V 



(S v -n u r]T), K 



A12 



N(0) 
W 2 
N(0)h 2 vl 
6W 2 



8n 2 T 2 
7C(3) ' 

(36) 



W 2 



where the energy scale W ~ ttT c is introduced for brevity 
and v v are the Fermi velocities in two bands which for 
simplicity are assumed isotropic. We, in fact, confirmed 
Eqs. (|36[) of Zhitomirsky and Dao employing different 
methods (except our b v is by a factor of 2 larger than 
that of Ref.Q). It is worth noting that the microscopi- 
cally derived a„ are not proportional to r as in the stan- 
dard one-band GL unless one of the parameters S v is 
zero; given the condition Q33p this may happen only if 
Ai 2 = 0. This feature of the two-band GL is sometimes 
overlooked i 21 ' 22 

As stressed in Ref.@, the term i^^jnA^ | 2 with order 
parameters gradients is the only possible in the GL en- 
ergy, although the symmetry may allow for other combi- 
nations of gradients. 

The coefficients entering the GL Eqs. (jTDJ) and pT|) are: 



N(0) 2 C 



K = 



V 



h 2 v 2 N(0) 2 
6W 2 r] 



Pu = 



N(0) 2 DS„ 

vw 2 x 2 2 



where 



v 2 



niS 2 v 2 + n 2 S\v\ 
has the dimension of a squared velocity and 

C = niS 2 + n 2 Si, D^niS 2 +n 2 S 2 



(37) 



(38) 



(39) 



are constants. 

Hence, we can express the length scale (fT8|) of the space 
variation of both Ai and A 2 in the GL domain in terms 
of microscopic parameters: 



2„~,2 



h 2 v 



2W 2 Ct 



(40) 



The upper critical field follows: H c2 = <Pq/2tt^ 2 . The 
one-band limit is obtained by setting m = 1, n 2 = 
so that C = S 2 and v 2 = S 2 v 2 /3 that yields £ 2 = 
7C(3)?iV/487r 2 T 2 T as it should. 

Variation of the free energy F with respect to the vec- 
tor potential A gives the current density. Following the 
standard procedure we obtain for the penetration depth 
of a weak magnetic field: 



1 

A 2 



32tt j 



E 



2„~,2 



167rCW(0)e 2 5 



(41) 



In the one-band limit this yields the correct result: A 
(167re 2 V(0)?; 2 /3c 2 )T. 
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A straightforward calculation yields the equilibrium 
zero-field free energy: 

F ^-N(0)W 2 ^t 2 . (42) 

The thermodynamic field H c follows: H 2 /8tt = —Fq. 
One can show that the relative specific heat jump at T c 
differs from the one-band value 12/7£(3) = f.43 by a 
factor C 2 /D < lM 

One can now form the dimensionless GL parameter, 



e 8irN(0)e 2 h 2 v 4 ' v ' 

and verify the standard relation H c ij H c \/2 = k. 

Finally, the equilibrium energy is evaluated by substi- 
tuting the solutions of the GL equations to the functional 
©: 

^=§J dV { b2 - 1 2^ 4 }- (44) 

where b = B/H c y2 is the dimensionless field. Thus, the 
theory of a two-band superconductor near T c is mapped 
onto the standard one-order parameter GL scheme. 

In particular, this mapping means that the GL prob- 
lem of the interface energy between normal and super- 
conducting phases has the same solution, i.e., n = 1/V2 
separates type-I and type-II superconductors. This has 
been demonstrated in Ref.0 by solving numerically the 
nonlinear system of GL equations (01,© without dis- 
carding terms 0(t 2 ) employed here. 

A. Remark on boundary conditions 

The solution (IT5|) for the two gap functions of the GL 
Eqs. (fTUj) and (fTTj) holds indeed provided the boundary 
conditions fo r A 2 a re the same as for Ai multiplied by 
the factor \fPifh- This is clearly the case for the ID 
problem of the S-N interface energy. The same is true 
for the problem of the single vortex structure: both A's 
are zero at the vortex center and approach A„ ; o with the 
correct ratio at infinity. 

However, for, e.g., proximity situations with a two- 
band superconductor on one side of the contact with a 
normal metal, the condition on the superconducting side 
far from the boundary is satisfied, whereas the question 
of boundary conditions at the boundary remains open. 
In this case, one cannot claim that A(a;)'s are propor- 
tional. Nevertheless, as is seen from Eqs. (fT0|) and (fTTj) . 
the length scale £ = yjK/ar is still the same for both 
order parameters. 

IV. DISCUSSION 

Two-band GL equation have been used in a number 
of publications where the coefficients in the GL energy 



functional a v ,b v , K v and 7 were varied and possible con- 
sequences were discussed. Moreover, different £'s and 
even A's were assigned to the two bands along with two 
different k's. This led to speculations that situations may 
exist where one of the bands behaves as a type-II super- 
conductor with Ki > l/y/2, while the other may have 
K2 < 1/V2 and behave near T c as the type-I; the su- 
perconductivity in such situations was called "type-1.5". 
MgB 2 has been suggested as such an example, see. e.g., 
and references therein. 

The present work argues that such situations do not 
exist. The point is that the GL equations are derived 
from the microscopic theory within certain approxima- 
tions that lead to the free energy near T c being propor- 
tional to (I — T/T c ) 2 and the order parameter (or param- 
eters) varying as (1 — T /Tc) 1 / 2 . Formally, the nonlinear 
system of GL equations (J4j) , d5j) for two-band materials 
can be solved with any accuracy. However, physically 
there is no point in going to accuracy higher than that 
of equations themselves; whatever results obtained along 
these lines will be absolutely unreliable. To get a near-T c 
description more accurate than GL, one should go back 
to microscopic theory that generates many extra terms 
in the free energy ex pan sion even for the one-band sit- 
uation, see, e.g., Ref.l24j. so that the multi-band gener- 
alization of such an approach is unlikely to produce a 
useful theory. It is demonstrated on a one-dimensional 
problem of Ref.@ and is shown for a general case in this 
paper that within the GL accuracy both order parameters 
of a two-band superconductor vary on the same length 
scale £ of Eq. (fT5| contrary to requirements of "1.5 type 
superconductivity" . 

We note that this conclusion holds for the "GL do- 
main" defined as the temperature interval near T c where 
the GL expansion can be justified. We do not specify 
this domain explicitly because its size may vary from 
one case to another. E.g., it is argued in Ref.@ that 
for two dirty bands (with no inter-band scattering) of 
MgB2, the domain of GL applicability shrinks practically 
to zero. However, whatever this size is, within this do- 
main the two order parameters vary on the same length 
scale. Therefore, attempts to employ the GL functionals 
- on hand - and to assume different length scales - on the 
other - cannot be justified. 

Moreover, we show that - within the GL accuracy - the 
two GL equations for the two-band case are reduced to 
a single equation for the normalized order parameter; in 
other words, the two-band superconductor is described 
by a single complex order parameter. This excludes pos- 
sibilities of having "fractional vortices" with exotic prop- 
erties such as those discussed in Refs.H^, [26l 

Microscopically, our results were derived within a 
weak coupling theory of clean superconductors. We 
believe, however, our conclusions go beyond that. For 
our results to hold it is crucial that due to the finite 
interband Josephson coupling 7, the coefficients a v in the 
GL energy remain finite at T c . Once this is guaranteed 
our qualitative conclusions remain unchanged, even if 
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assumptions of the weak coupling, no-scattering, and 
isotropy do not apply. 

Discussions with R. Fernandes, J. Geyer, J. Clem, 
R. Prozorov, M. Das, M. Milosevic, M. Zhitomirsky, 
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